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Abstract. We present a new approach to obtaining the scaling behavior of the 
entanglement entropy in fractional quantum Hall states from finite-size wavefunctions. 
By employing the torus geometry and the fact that the torus aspect ratio can be readily 
varied, we can extract the entanglement entropy of a spatial block as a continuous 
function of the block boundary length. This approach allows us to extract the 
topological entanglement entropy with an accuracy superior to what is possible on 
the spherical or disc geometry, where no natural continuously variable parameter is 
available. 

Other than the topological information, the study of entanglement scaling is also 
useful as an indicator of the difficulty posed by fractional quantum Hall states for 
various numerical techniques. 
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1. Introduction 

Describing condensed matter phases using entanglement quantifiers from quantum 
information theory is a rapidly growing interdisciplinary topic [1, 2]. Of special interest 
are the relatively rare cases where entanglement can provide information not readily 
captured by conventional quantities such as correlation functions. This is the situation 
for systems possessing topological order [3], where entanglement considerations have 
proven useful [4, 5, 6]. In particular this has led to insight into the structure of fractional 
quantum Hall (FQH) states [7, 8, 6, 9, 10, 11, 12, 13, 14]. In this Article we focus on 
entanglement in this most realistic class of topologically ordered states. 

A prominent measure of entanglement is the von Neumann entropy of entanglement, 
Sa, measuring the entanglement between a block (A) and the rest (B) of a many-particle 
system in a pure state. The entanglement entropy Sa = — tr [pa In pa] is defined in terms 
of the reduced density matrix, pa — tr B p, obtained by tracing out B degrees of freedom 
from the system density matrix p = \ip)(tfj\, with denoting a ground state wave 
function. In one dimension the scaling behavior of the block entanglement entropy is 
well understood, see, e.g., Ref. [2]. In two dimensions (2D), no such generic classification 
exists. However, for topologically ordered states in two dimensions, the entanglement 
entropy contains topological information about the state: Sa scales as 

S A = aL-nj + 0(l/L), (1) 

where L is the block boundary length, 7 characterizes the topological field theory 
describing the state [4, 5] , while n counts the number of disconnected components of the 
boundary. The value of 7 is related to the "quantum dimensions" of the quasiparticle 
types of the theory, 7 = InP, where T> is the total quantum dimension. For Laughlin 
states at filling v = 1/m, 7 = | lnm. For more intricate FQH states, some examples of 
7 values are provided in Refs. [4, 15, 8]. If 7 can be determined accurately, its value can 
in principle be used to determine whether a topological phase belongs to the universality 
class of a given topological field theory. 

A numerical determination of 7 and a requires information about Sa for a number 
of different boundary lengths, L. In Ref. [7] such information was obtained from finite- 
size FQH wavefunctions by approximating the spatial partitioning by partitioning of 
the discrete set of Landau level orbitals. Refs. [7, 8, 10] used spherical geometries, and 
explored ways of extrapolating entanglement information from such geometries to the 
thermodynamic limit. Ref. [12] used disk geometries to calculate 7 for bosonic FQH 
wavefunctions. The accuracy in the determination of 7 from finite-size wavefunctions 
on these geometries remains disappointing (10%-30% for the simplest Laughlin states). 
Improved methods for calculating 7 are thus sorely needed. 

In this work, we report a significant advance in this direction, through the use of 
the torus geometry and the fact that the aspect ratio (circumference) of the torus can be 
varied continuously without drastically altering the torus setup or symmetry. Varying 
the circumference changes the length of the boundary between A and B. No natural 
analogous continuous parameter exists in the other geometries, so that in those cases 
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each system size and bipartition provides only one point in parameter space. Exploiting 
the continuous parameter, we present a procedure that leads to an accuracy in 7 down 
to a few percent. Our analysis also provides a visual and physical indication of the 
reliability of the extracted 7 value. 

Previously, Ref. [16] reported a naive modification of the sphere algorithm of Ref. [7] 
to the torus geometry with fixed aspect ratio. We will show why this analysis was 
inappropriate and based on an extrapolation procedure that is not meaningful for the 
torus geometry. 

In addition to the topological content of the subleading term in (1), the dominant 
linear term itself is also of some importance. The rate of entanglement growth, a, 
indicates how challenging the state is to simulate on a classical computer, through a 
one-dimensional algorithm like DMRG [17, 18], or through recently-proposed true two- 
dimensional algorithms like PEPS [19] or MERA [20]. DMRG has been used to simulate 
FQH states [21, 22, 23, 24], and these states pose a future challenge for two-dimensional 
algorithms currently under development. 

The calculation of the topological entanglement entropy 7 is of significant current 
interest, not only for FQH states but also for various other topologically ordered states. 
For the zero-temperature Kitaev model, it is relatively easy to calculate 7; so the 
concept has been used in exploring issues such as temperature effects and quantum phase 
transitions [25, 26, 27, 28]. For quantum dimer models and related states, considerations 
of entanglement scaling are more intricate [29, 30], of difficulty comparable to FQH 
states. More generally, entanglement scaling in 2D states of all kinds has become the 
focus of intense study at present [31, 2, 32]. 

In Section 2 we show how the torus geometry allows us to map the interacting 
Landau level (LL) problem onto a one-dimensional lattice model, appropriate for 
studying bipartite entanglement. In Section 3 we outline the general behavior 
of entanglement entropy on the torus geometry and deal with the issue of torus 
degeneracies. In Section 4 we present our main results, including an analysis leading 
to the determination of 7. The concluding Section 5 connects to the existing literature 
and discusses implications of our results. 

2. Torus setup — geometry, partitioning 

We study an A-electron system on a torus (see Fig. 1) with periods Li,L 2 in the x- 
and ^-directions, satisfying LiL 2 = 2ttN s in units of the magnetic length. The integer 
N s = N/v is the number of magnetic flux quanta. In the Landau gauge, A = BySt, a 
basis of single particle states in the lowest Landau level can be taken as 

m 

with j = 0, 1, N s — 1. The states ipj are centered along the lines y = —2-nj/Li. Thus 
the ^/-position is given by the x-momentum j. 
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Figure 1. Geometry of the torus and bipartitioning. The lowest Landau level 
is spanned by orbitals which in the Landau gauge are centered along the circles shown. 
On the right, we represent the torus as a rectangular region with periodic boundary 
conditions in both directions. The dimensions of this rectangle (Zi, L 2 ) are the two 
circumferences of the torus. The example shown here has N s = 12 orbitals with Ia = 4 
orbitals in the A block. 



A generic translation-invariant two-body interaction Hamiltonian, acting within a 
Landau level, can be written as 



H 



n fc>|m| 



(3) 



where creates an electron in the state ip m and Vk m is the amplitude for two particles 
to hop symmetrically from separation k + m to k — m [33]. Hence, the problem of 
interacting electrons in a Landau level maps onto a one-dimensional, center-of-mass 
conserving, lattice model with lattice constant In jL\. This provides a natural setting 
for defining entanglement, by bipartitioning the system into blocks A and B, which 
consist respectively of I a consecutive orbitals and the remaining Ib = N s — I a orbitals 
(Fig. 1). Since the orbitals are localized in the direction of the lattice, this is a reasonable 
approximation to spatial partitioning, as on the sphere [7, 8, 6, 10, 11]. 

Because this partitioning implies two disjoint edges between the blocks, each of 
length Li, the entanglement entropy should satisfy the following specific form of (1): 

S A (L 1 ) = 2aL 1 - 2 7 + 0(l/L 1 ). (4) 

Thus our setup should yield a linear scaling form of the entropy with the L\ = intercept 
at -27. 

In this work, we obtain ground states of (3), in the orbital basis (2) using the 
Lanczos algorithm for numerical diagonalization. We study bipartite entanglement in 
these ground states. Apart from diagonalizing the Coulomb problem we also consider 
pseudopotential interactions [34, 35] which have the Laughlin states [36] as exact ground 
states. The largest Hilbert space sizes considered are 208'267'320 for 39 orbitals at 
v = 1/3, 19'692'535 for 45 orbitals at v = 1/5 and 66'284'555 for 35 orbitals at v = 2/5. 
The simulations are however currently limited by the size of the reduced density matrices 
to be calculated and fully diagonalized. 
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Fractional quantum Hall states have degenerate ground states on the torus 
geometry. It is convenient to label the ground states by their corresponding thin torus 
(or Tao-Thouless, TT) patterns [33, 37, 38, 39, 40]. For example, for v = 1/3 there are 
three degenerate states, which correspond to the TT configurations 

100100100 



100100100 
010010010 
001001001 



100100100100100100 
010010010010010010 



poiooiooiooiooiooi 



010010010 

001001001 , (5) 



A 

for N s = 36. Here the positions of l's indicate the positions (or, equivalently, the 
transverse momenta) of filled single particle states. An equal partitioning {I a — Ib — 
N s /2) is illustrated. 

In general, abelian FQH states at v = p/q have q degenerate ground states, related 
to each other through translation and corresponding to q thin-torus patterns, each 
composed of unit cells with p electrons on q sites. These states are ground states 
for generic (two-body) interactions as L\ — > [33]. For non-abelian states there is 
an enhanced degeneracy and the corresponding thin torus patterns are not simply 
translations of each other. 

The thin torus states are unentangled product states, in the orbital basis. As 
Li is increased from zero, fluctuations on top of the TT states will make the states 
entangled. A crucial property of the FQH states is that their bulk versions are, for 
appropriate interactions, adiabatically connected to their respective TT states and the 
gap is finite for all L\ [33, 37]. This allows us to probe the response of these states 
as the geometry is deformed. Such deformations have also been considered to extract 
properties such as the Hall viscosity [41, 42], to put consistency conditions of FQH states 
[43], to find instabilities to competing states (see eg [44, 45]) and to deform the torus 
to the solvable thin torus limit as discussed above. All the three fractions studied in 
this paper [y = 1/3, 1/5, and 2/5) are, for the pseudopotential as well as the Coulomb 
interactions, continuously connected to their TT states. 



3. Degeneracy averaging, Area law at constant L\ 

For any finite L\ the charge density modulations of the TT pattern will prevail to some 
extent, leading to different entanglement in the q degenerate ground states. We illustrate 
this in Fig. 2(a) where we plot the entanglement entropy Sa{Ia) as a function of Ia in 
the three degenerate v = 1/3 Laughlin wave functions. For each I a two out of the three 
entanglement entropies are equal, while the third one is different, as can be inferred 
from examining the partitionings shown in Equation (5). Two of the TT patterns have 
1-0 and 0-1 cuts at the block boundaries, while the third has only 0-0 cuts. Since the 
microscopic environment at the two boundaries determines the entanglement spectrum 
[13] and hence the entanglement entropy, this implies that two of the entanglement 
entropies are equal at each I a- 
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Figure 2. Degeneracy averaging and I a dependence, (a) Entanglement entropy 
in different degenerate sectors for the v = 1/3 Laughlin state, and their arithmetic 
average, (b) difference between degeneracy-averaged Sa and largest individual-sector 
Sa, as a function of L±. They differ significantly only for intermediate L±. (c) 
Degeneracy-averaged Sa versus I a, for different L\ values. 

The Sa{Ia) each have prominent oscillatory behavior. However, we find that the 
arithmetic mean of the three individual entanglement entropies is remarkably free 
of oscillations. We will thus base all our following discussions on the degeneracy- 
averaged entropy J. Ultimately this averaging will become unimportant for very large 
Li, as Fig. 2(b) shows that the difference AS between the maximum and the mean 
entanglement entropies in the q sectors vanishes rapidly at large L\ (starting to decrease 
after L\ ~ 10 in the Laughlin v = 1/3 case shown). 

In Fig. 2(c) we study the I a and N s dependence for a given L x . At constant L x we 
expect the entropy to saturate once I a is large enough, since the block boundary length 
(2Li) is held constant. This is indeed what is found numerically in Fig. 2(c). The length 
scale controlling the saturation is the (real space) correlation length £ r in the y direction 
of the incompressible FQH liquid. The correlation length £ Q measured in number of 
orbitals is expected to scale as ~ £ r x Li/2n. The saturation of the entanglement 
entropy for large I a is in complete analogy to the area law for one-dimensional gapped 
systems [46, 47]. It is this saturation value Sa{L\) of the entanglement entropy obtained 

\ In principle the averaging could be done in other ways, for example, one could average over the 
density matrices p — > ^2^=1 P^V^gs or reduced density matrices, and then compute the entanglement 
entropy from this averaged matrix. We do not however pursue such alternate averaging procedures in 
the present work. 
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Figure 3. v = 1/3 Laughlin state: Entanglement entropy, (a) Sa and (b) its 

derivative dSU/dLi for the Laughlin state at v = 1/3 as a function of L\. From the 
plateau behavior in dSU/dii for L\ > 12 we infer a s» 0.153(2). 

for I a 3> £ that will be analyzed in the following. To avoid the finite size effect as far 
as possible we consider I a = N s /2 for N s even, and I a = (As — l)/2 for N s odd, in the 
rest of this article. From Fig. 2(c) we can again infer that the averaged Sa indeed has 
a much smoother dependence on the block size compared to the Sa for the individual 
degenerate states. 

4. Accessing the scaling regime 

In this Section we provide our main numerical results and discuss how a and 7 can be 
extracted by continuously varying L\. As mentioned above, from now on Sa refers to the 
equal-partitioning entanglement entropy (I a = N s /2 for N s even, and I a = (N s — l)/2 
for N s odd). 

v = 1/3 Laughlin state - Fig. 3 shows the behavior of Sa{L\) (a) and its derivative 
oISa/oILi (b), for the Laughlin state at fraction v = 1/3, arguably the most prominent 
and also the simplest FQH state. In this figure and subsequently, we use a five-point 
formula to numerically obtain derivatives. The entanglement entropy Sa(Li) remains 
minuscule until L\ ~ 3, and then gradually changes to the expected linear increase 
behavior which is reached around L\ ~ 7. There are oscillations on top of the linear 
behavior, which are more prominent in the derivative plot (b). The oscillations can be 
interpreted as an interplay between the finite circumference along the x direction (finite 
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Figure 4. v = 1/3 Coulomb ground state: Entanglement entropy, (a) S A and 

(b) its derivative dSU/dLi for the Coulomb ground state at v = 1/3 as a function of 
L\. One curve for the Laughlin state is also shown for comparison (dashed line). In the 
inset (c) the difference in Sa (Li) between the Coulomb and the Laughlin is displayed 
(for N s = 30,36). 

Li) and the interparticle distance. The oscillations die off as a function of Li, so that if 
N s is large enough one can get the scaling form at large L\. 

For small L\ the finite size convergence is essentially perfect. At larger L 1; the 
Sa{Li) and dS A /dL\ curves show stronger dependence on N s . The ^-dependence 
shows up first for the smallest system sizes and at increasing L\ for progressively larger 
system sizes. This reflects the fact that, for any finite-size system, at very large L\ the 
edges of A get too close (small L 2 ) and cannot be thought of as independent [13]. In 
particular, once L\ exceeds some value we enter the "dual thin torus" or "thick torus" 
limit [40, 48, 43], and the entanglement entropy levels off to some saturation value. 
Corresponding to the saturation of Sa(Li), the derivative oISa/oILi drops off to zero 
after some L\. 

Thus, the scaling form of (4) is valid only in a window of Li, after the oscillations 
have subsided but before S A {L\) saturates, or shows other precursor finite size effects. 
This plateau region can be seen clearly in the oISa/oILi curve for the v = 1/3 Laughlin 
data for L\ > 12. The finite size convergence of the data also provides a clear signal 
showing whether the bulk scaling regime is reached or whether (geometrical) finite size 
effects are still significant in the numerically accessible L\ regime. 
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Figure 5. v = 1/5 Laughlin vs Coulomb: Entanglement entropy (a) and its 
derivative (b) for the v = 1/5 Laughlin wave functions. Entanglement entropy (c) and 
its derivative (d) for the v — 1/5 Coulomb ground state. 



v = 1/3 Coulomb ground state - Fig. 4 plots (a) Sa(Li) and (b) oISa/oILi for the 
Coulomb ground states at v — 1/3. While this state has somewhat more severe finite- 
size effects and oscillatory behaviors compared to the model Laughlin state of Fig. 3, 
we note that the scaling form of their entanglement entropies are very similar. To 
further highlight this fact we plot the difference between the entanglement entropies as 
a function of Li in the inset of Fig. 4(c). The similarity of entanglement entropies of the 
two states is not unexpected from the perspective that the states have a large overlap 
for all Li [45], but it is nevertheless interesting considering that a more "generic" state 
such as the Coulomb state is expected to have larger entanglement. One could thus 
have expected the Coulomb state to have a larger a, as defined in Eq. 1, but Fig. 4(b) 
suggests a very similar a ~ 0.15(1). 

v = 1/5 - Fig. 5 shows the Sa(Li) and oISa/oILi behaviors at v — 1/5, for both the 
Laughlin (a) & (b) and the Coulomb ground state (c) & (d). As expected, the finite-size 
oscillations are much more severe in these states. This is expected as the interparticle 
distance is larger; thus larger systems should be required to reach the scaling regime. 
Moreover, the proximity to the Wigner crystal phase makes the Coulomb ground state 
deviate more substantially from the Laughlin state than is the case at v = 1/3 [49, 45]. 
While we are able to get an almost ^-converged S(Lx) curve up to L\ ~ 18 for the 
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Figure 6. Li-local extraction of 7. The intercept of local linear approximations 
to the Sa{La) curves, i.e., Sa{L\) — L\ x dS/dLi, plotted as a function of L\. In the 
scaling regime, this quantity should give —27. The symbols for v = 1/3 (y — 1/5) 
are the same as in Fig. 4 (Fig. 5). Theoretically expected —27 values are shown as 
dashed horizontal lines. In panel (c) the solid line through the largest size data is the 
fit obtained using Equation (6). 



Laughlin state [leading to a rough estimate of a ~ 0.17(2)], the finite size effects in the 
Coulomb ground state are so severe that no meaningful extraction of a is possible with 
current system sizes. 

Extraction of the topological entanglement entropy 7 - In Fig. 6 we show calculations 
of 7 for the Laughlin state at v — 1/3 (a) and v — 1/5 (c) as well as for the Coulomb 
ground state at the same fractions (b) & (d). Evaluating the L 1 derivative using a 
centered 5-point formula, we plot SA{Li)—LixdS/dLi as a function of L\. This quantity 
is the intercept of a linear approximation made to the Sa{L\) curve locally at each L\. It 
should take the value —27 in the scaling region, see Equation (4). Not surprisingly, the 
intercept oscillates at intermediate Li, has a plateau in the "scaling window" described 
above, and then moves off to a large positive value when Li is yet larger entering the thick 
torus regime. The plateau region value gives us the best estimate for the topological 
entanglement entropy. A significant advantage of our analysis is that, by examining 
the oISa/oILi curve (and its N s dependence), we can identify the correct window of L\ 
values to use. 

In Fig. 6(a) the v = 1/3 Laughlin shows such a clear plateau region. The 
plateau region value around L\ ~ 18 gives us the best estimate for the topological 
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Figure 7. v = 2/5 Coulomb ground state, (a) Sa and (b) its derivative dSU/dZa 
for the Coulomb ground state at v = 2/5 as a function of L\. (c) Estimation of —27 
through a plot of the Li-local intercept against L\. 



entanglement entropy 7 ~ 0.565(5), to be compared to the theoretical expectation 
7 = ln(3)/2 pa 0.5493. The difference amounts to only 3 percent in this ideal case. 

For the u = 1/5 Laughlin [Fig. 6(c)], the finite-size issues are significantly larger, 
and the oscillations have not yet damped out at accessible sizes. However, one can take 
the average of the oscillating values to get a reasonable estimate of 7. We use a simple 
damped oscillation fitting ansatz of the form: 

f(Li) = —27 + a x exp[— bLi] x sin(cLi — d), (6) 

and fit the N s = 45 curve for L\ > 9, yielding an estimate of 7 ~ 0.81. This value again 
compares very favorably to the theoretical expectation 7 = ln(5)/2 pa 0.8047. 

At each of these fractions, the finite-size convergence is worse for the Coulomb 
ground state compared to the Laughlin ground state. While for the v = 1/3 Coulomb 
[Fig. 6(b)] a fitting analysis along the lines of the v — 1/5 Laughlin still provides a 
reasonable 7 estimate: 7 ~ 0.60, the v = 1/5 Coulomb state [Fig. 6(d)] clearly does not 
allow a meaningful 7 extraction from system sizes presently reachable through numerical 
exact diagonalization. 

v — 2/5 Coulomb interaction - Finally, in Figure 7 we consider the Coulomb ground 
state at filling v — 2/5, whose 7 value has not been studied numerically so far. This state 
is, for all L\, well described by the torus version [50] of the Jain [51] (or, equivalently, the 
hierarchy [34]) state. The finite-size effects are somewhat less severe than the v = 1/5 
Coulomb case [Figs. 5(c,d) and 6(d)]. One obtains an entanglement growth rate of 
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a ~ 0.188(16). While the ^-convergence is not good enough for a precise determination 
of 7 (expected to be |ln5), examination of the largest two available sizes suggests that 
two or three additional sizes may be enough to provide an estimate at the ~10% accuracy 
level. 

5. Discussion 

In this article we have shown how continuous geometric deformations of the torus can be 
employed to explore the scaling form of the entanglement entropy. This has allowed us to 
propose a method for determining the topological part, 7, from finite-size wavef unctions, 
to greater precision compared to earlier analyses which did not utilize any continuous 
parameter. Our analysis indicates that current state-of-the-art system sizes are enough 
to obtain reliable 7 calculations for the simplest fractional quantum Hall states (Laughlin 
states), but that more intricate states would require larger sizes than currently accessible, 
in order to reach the scaling limit. Our procedure provides a clear method for identifying 
whether the scaling window has been reached or not. 

There has been an earlier report of entanglement entropy and 7 calculations on 
the torus [16], using fixed aspect ratio, Li/L 2 = 1. Ref. [16] performed N s — > 00 
extrapolations at fixed I a, and expected the extrapolated values to scale as c\\/Ia — 27. 
We illustrate such a fixed I a extrapolation in Fig. 8. The extrapolation does not lead 
to a physically meaningful limit because the boundary lengths diverge and the two 
boundaries get infinitesimally close to each other, in the N s — > 00 limit. 



BmB 





Figure 8. Extrapolation at fixed aspect ratio and fixed I a- Fixed I a implies 
that the area covered by the region A is constant. Such extrapolation does not lead to 
a well-defined limit as the limiting case is one with infinitesimal thickness and infinite 
boundary length. 

The extrapolation of Ref. [16] arose from an incorrect adaptation of the procedure 
of Ref. [7] which was designed for the sphere geometry. For the spherical case, the 
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fixed-Z^ extrapolation takes one to the well-defined limit of an infinite disk with circular 
A region, where the boundary length indeed scales as ~ v^a- Also for the disk [12], 
the fixed-^A extrapolation to N — > oo is well-defined as the limiting situation is again 
an infinite disk with constant density. (The density is not constant in finite-size disk 
simulations.) However, for the torus with unit aspect ratio, the fixed-Z^ extrapolation 
to N — > oo has a pathological limit for the shape of the A block, and the limiting 
entanglement entropy has no reason to scale as ~ \[Ta- A plot of the extrapolated Sa 
versus \[Ta thus has no obvious connection to the entropy scaling as a function of the 
boundary length, or to the definition of 7 as formulated in Refs. [4, 5]. 

The idea of obtaining entanglement entropy scaling through varying discrete torus 
circumferences has been employed in Ref. [29] for the dimer model on the triangular 
lattice. The details are quite different from the FQH case. It is possible that some of the 
ideas developed here for the FQH context might be transferred fruitfully to numerical 
work on dimer models or other lattice models. Ref. [32] has studied entanglement of 
integer quantum Halls states on a torus, between true spatial partitions rather than 
orbital partitions. The orbital partitioning entanglement is zero for integer quantum 
Hall states because they are product states in the orbital basis. 

Since Ref. [7, 8] reported entanglements of the same Laughlin states on a different 
geometry, it is interesting to compare the magnitudes of the entanglement entropy. The 
data tabulated in Ref. [8] for the v — 1/3 state yields a value of a ~ 0.15, which is 
close to that obtained from the torus entanglement data reported in this work. While 
this is not unexpected, it has several instructive implications. First, it can be regarded 
as additional evidence that orbital partitioning entanglement is a good approximation 
to spatial partitioning entanglement. Second, it shows that the entanglement entropy 
contributions from the two edges simply add for a block with two edges compared to 
one edge. 

In addition, the fact that the sphere and the torus have the same "entanglement 
entropy density" per unit boundary length, allows us to compare the difficulty of DMRG 
simulations on spherical and toroidal geometries based on considerations of linear sizes 
alone. Conventional wisdom might be that the torus is significantly more difficult to 
treat using DMRG, because it is a system with periodic boundary conditions while the 
sphere is more analogous to a lattice system with open boundary conditions. However, 
the torus has the same block boundary length everywhere (2L\ = 2\/2ttN s for the 
hardest case of unit aspect ratio), while the block boundary on a sphere varies. For 
DMRG on a spherical geometry, the dominant reduced entropy contributions comes 
from the equator region, where the block boundary is L = TTy/2N s . The torus boundary 
is only a factor 2/^/7? rs 1.13 larger than the sphere case, rather than being twice as 
large. Moreover, from the two edges one benefits twice the (negative) contribution from 
the topological part of the entanglement entropy (in contrast to a single contribution on 
the sphere). We therefore infer that DMRG on the torus geometry is not as drastically 
more difficult compared to the sphere would be suggested by the argument of 

two block boundaries versus one. 
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The entanglement entropy S in a bipartition of a quantum state puts a lower bound 
on the number of states m to be kept in an accurate DMRG simulation through the 
relation m ~ exp(S'). The computational complexity of DMRG being polynomial in m, 
the scaling of S with system geometry is thus of primary importance. In the simplest 
cases, e.g., one dimensional gapped quantum systems with local interactions, the entropy 
S does not depend on the block length, enabling DMRG simulations for basically infinite 
systems at constant m. Torus FQH simulations at constant L\ [22] also belong to this 
tractable class. If one is however interested in describing true bulk FQH systems at 
fixed aspect ratio (Li/L 2 = const), then entropy S will scale linearly with Li oc \^N\ ! . 
This translates into m oc exp [const x \/N s \, i.e., accurate DMRG simulations for bulk 
FQH states scale exponentially in the physical width, similar to 2D lattice models [18]. 

The present work opens up a number of directions deserving exploration. Our 
analysis provides a way to decide on whether or not available wavefunction sizes provide 
access to the entanglement scaling regime. Thus, as bigger wavefunctions become 
available, our analysis can be applied directly to obtain better calculations of the 
topological entanglement entropy 7. Eventually, this type of calculation could become 
a standard tool for diagnosing unknown or poorly-understood FQH states. Another 
obvious direction is the study of entanglement scaling through geometric deformations at 
other fractions and more complicated states. It may also be interesting to try to devise 
continuous geometric tuning parameters for other geometries, e.g., one may consider 
ellipsoidal geometries as deformations of the sphere, although setting up the Landau 
level problem in such geometries is not straightforward or convenient. It is possible 
that a combination of various deformation considerations might lead to further refined 
procedures for estimating entanglement scaling and 7. 
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